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Abstract 

We show that a parametrically driven cubic-quintic complex Ginzburg-Landau equation exhibits a hysteretic nonequilib- 
rium Ising-Bloch transition for large enough quintic nonlinearity. These results help to understand the recent experimental 
observation of this pheomenon [A. Esteban-Martm et al, Phys. Rev. Lett. 94, 223903 (2005)]. 



PACS numbers: 42.65.Sf, 47.54.+r,42.65.Hw 

Spatially extended bistable systems with broken phase 
invariance display defects in the form of interfaces, so- 
called domail walls (DWs). A paradigm for the study 
of DWs is the parametrically driven complex Ginzbur 
Landau equation, which can be written in the form 



[ilf - 



d t A = 1 A* + {p + iv)A + {l + ia) d 2 x A - (1 + i[3) \A\ 2 A, 

where 7 is the parametric pump, /i accounts for linear 
gain or loss, depending on its sign, v is a detuning, a 
is the diffraction coefficient, and (3 is the nonlinear dis- 
persion coefficient. In writing Eq. the spatial co- 
ordinate and the field amplitude have been normalized 
to the square root of the diffusion coefficient and of the 
saturation coefficient, respectively. This equation rep- 
resents a universal description of parametrically excited 
waves as well as of the close to threshold dynamics 
of self-oscillatory systems externally forced at the second 
harmonic of the natural oscillations frequency . 

In Eq. Q the phase invariance of the field A (x, t) is 
broken because of the presence of the parametric term 
-fA*, i.e, Eq. Q shows the discrete symmetry A «-> —A. 
This makes possible the existence of domain walls (DWs) 
that connect spatial regions where the field passes from, 
e.g., the homogeneous solution Aq to the equivalent so- 
lution — Aq. There are two types of DWs, namely Ising 
and Bloch walls, which differ in the way the field crosses 
the complex zero at the DW core: In the Ising wall, both 
the real and the imaginary parts of the field become null, 
whilst in the Bloch wall the real and the imaginary parts 
become null at different spatial points. Then, in terms of 
the field intensity, |^4| 2 , an Ising wall is dark at its center 
whilst it is grey in the case of a Bloch wall. But the most 
striking difference between Ising and Bloch walls lies in 
their different dynamic behaviour: When nonvariational 
terms are present (in Eq. Q this means that v, a or 
(5 be different from zero) Bloch walls move whilst Ising 
walls remain at rest. 

Coullet et al. 0] have discussed the above in detail 
for fi > and have shown how Ising walls bifurcate into 
Bloch walls through the so-called nonequilibrium Ising- 
Bloch transition (NIBT), that takes its name from the 



equilibrium Ising-Bloch transition ocurring in ferromag- 
nets 0. Subsequently it was also shown that the NIBT 
occurs in Eq. for negative ijl 0, in which case DWs 
connect not only equivalent homogeneous solutions but 
also spatially modulated solutions. 

There are a few experimental observations of this phe- 
nomenon. As far as we know, it has been reported only in 
liquid crystals Q either subjected to rotating magnetic 
fields, 13,13 or to an alternate electrical voltage This 
last experiment constitutes a particularly clear observa- 
tion of the NIBT free from 2D effects, which complicate 
front dynamics through curvature effects. We must add 
our very recent observation of a hysteretic Ising-Bloch 
transition in a nonlinear optical cavity ||. This last ex- 
periment was carried out in a photorefractive oscillator 
in a degenerate four-wave mixing configuration 0, 0, Hcj 
and the cavity detuning played the role of the control 
parameter. We found that for small positive cavity de- 
tuning the system exhibits Ising walls. When detuning 
was increased, Ising walls bifurcated into Bloch walls at 
a cavity detuning value vib, and for v > v\b DWs were 
always of Bloch type. Then, when making a reverse de- 
tuning scan we found that Bloch walls existed up to a 
detunig value vbi where a new Ising-Bloch transition 
occurs. The interesting thing is that vbi < vib, and 
then there is a detuning domain, vbi < v < vib, where 
Ising and Bloch walls coexist. 

The origin of the hysteresys was experimentally found 
to lie in the existence of bistability in the homogeneous 
state of the system: The homogeneous state exhibits 
bistability within a certain cavity detuning range be- 
tween two spatially homogeneous states, say Ahom.i and 
Mxomfl, with lAhom^l 2 > |A horo)2 | . It occurs to happen 
that DWs connecting Ahom,i with — ^4hom,i are of Ising 
type whilst those connecting Ah om ,2 with — A^ om< 2 are 
of Bloch type. Eq. does not give any insight into 
this type of behaviour as the NIBT it exhibits is not hys- 
teretic nor its homogeneous solution exhibits bistability 
of the type we are describing. Here we try to put some 
light into this problem by considering a straightforward 
generalization of Eq. 0J. 
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Hysteretic Ising-Bloch transitions have been theoret- 
ically described recently in two very different contexts. 
On the one hand, it has been described in the anisotropic 
XY-spin system in an oscillatory magnetic field [JjJ . In 
that paper, Eq. (JIJ is studied for v — a — — plus 
an additional modulation term. As the system under 
study is variational, Bloch walls do not move and the 
hysteresys is found on the average oscillation period of 
the DW for certain parameter sets. On the other hand, 
hysteresis has also been found in the routes leading from 
standing fronts to a couple of counterpropagating ones 
in two discrete models (an array of Lorenz units and the 
FitzHugh-Nagumo model) in Ref.^3]. These previous 
theoretical results do not help to understand the exper- 
imental results in M that we have just resumed. As al- 
ready outlined in |8( we will show that the addition of 
a quintic nonlinearity in Eq. Q allows to understand 
qualitatively the experimental results. 

We start with the following natural generalization of 
Eq. Q 

d t A = j A* + ( A t + iv) A + (1 + id) d 2 A 

-(l + i0 3 )\A\ 2 A-(a 5 +i0 5 )\A\ i A, (2) 

i.e., we have added a fifth order nonlinearity with as and 
0s the quintic saturation and nonlinear dispersion coeffi- 
cientes, respectively. As we are interested in the minimal 
modification of Eq. JTJ that contains a hysteretic NIBT, 
in the following we shall concentrate on the special case 
as = 0, as the addition of the quintic nonlinear dispersion 
term is enough for our purposes as we show next. 

In Fig. 1 we represent the square modulus of the ho- 
mogeneous solution of Eq. © as a function of detun- 
ing v for 7 = 2, fj, = a = — 03 = 1 and the different 
values of 0s indicated in the figure inset. Notice that 
there are two regions in which the homogeneous solu- 
tion is multivalued: For negative detuning, where there 
is coexistence between two homogeneous solution values 
and the trivial solution, and also for positive detuning 
whenever /3s > /3g ~ 0.39428, where there is coexistence 
of three homogeneous solutions. We find that the lat- 
ter requires that the signs of the nonlinear dispersion 
coefficients, 03 and 0s are different, and will concen- 
trate on this case. Note that Eq. (0 holds the sym- 
metry (A,L>,a,P 3 ,f3 5 ) < — ► (A*, -v, -a, -03, -0s) and 
consequently the behaviour of the homogeneous solu- 
tion is the same for the parameter sets (y, 03,0s) and 

{-V, ~0 3 , -0s). 

The numerical integration of Eq. J3J shows that for 
negative detuning v, the system exhibits extended pat- 
terns and that the homogeneous solution can be stable 
for positive v. Then it is for v > that we can find DWs 
connecting homogeneous solutions and we concentrate in 
this case (Notice that for 5 = this is the parameter 
region, v > 0, where the NIBT was studied in Q] and 



We have carried out the numerical integration of Eq. 
(j2J for 7 = 2, /J, = a = 1 and different values of 03 and 
0S- We pass to comment first our results for 03 = —1 
and different values of 0s and v. 

In Fig. 2 we represent again the homogeneous steady 
state for the same parameters as in Fig. 1 (except for 
the values of 0s), and have marked the different patterns 
one can observe. For 5 — 0.38 < 0%, Fig. 2(a), the 
homogeneous solution is single-valued and two types of 
DWs are observed: Ising walls (IW in the figure), for 
detunings v < 1.24; and Bloch walls (BW) for v > 1.24. 
For v > 1.53 the homogeneous solution becomes modu- 
lationally unstable; Bloch walls connect patterns in this 
region. In Fig. 3 both the intensity and phase spatial 
profiles corresponding to an Ising wall, Fig. 3(a), and a 
Bloch wall, Fig. 3(b), are represented. Notice that the 
field intensity is null at the DW core only in the Ising 
wall, and that the phase jump is sharp (smooth) for the 
Ising (Bloch) wall. It is also interesting to notice the 
shoulder in the field intensity in the case of the Bloch 
wall (the shoulder appears on the back side of the wall 
with respect to the direction of movement). 

As the value of 0s is increased, new features appear. 
For 05 = 0.395 (i.e., slightly larger that 0$) the homo- 
geneous solution becomes multivalued, Fig. 2(b) (the 
dashed line indicates that the homogeneous solution is 
unstable). In this case Ising walls do not bifurcate di- 
rectly into moving Bloch walls, but start oscillating pe- 
riodically around a fixed position, that is, there is not a 
net displacement of the wall. Fig. 4 shows the intensity 
and phase profiles of the oscillating wall in three different 
instants of time and it can be appreciated how the wall 
pases from a clear Bloch character (smooth phase jump) 
when it is at the center of the oscillation, Fig. 4(b), to 
a clear Ising character (sharp phase jump) when it is at 
the extremes of the oscillation, Figs. 4(a) and (c). 

The behaviour just described appears when 0s — 0%. 
When 0s is further increased a new and remarkabe effect 
appears: There is a detuning range of coexistence be- 
tween the oscillating and the Bloch walls, see Fig. 2(c). 
This is better appreciated in Fig. 5 where the velocity of 
the walls is represented as a function of detuning for the 
same parameters as in Fig. 2. In Figs. 5(a) and 5(b) the 
behaviour of the velocity closely follows that of the stan- 
dard NIBT but in Fig. 5(c) the new phenomenon 
of the hysteretic NIBT is clearly appreciated. 

In the case we have described no coexistence between 
Ising and Bloch walls is observed, only between oscillat- 
ing and Bloch walls. But by decreasing the value of 03 
from —1 to —1.5 we can observe this coexistence. In Fig. 
6 we represent again the wall velocity as a function of de- 
tuning for the same values as in Fig. 2 except 03 = —1.5 
and 0s = 0.6. The behaviour is similar to that described 
above but now there appears a wide domain of coexis- 
tence between Bloch walls and both oscillating and Ising 
walls. This is a much clearer hysteretic NIBT. 
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These results show that the addition of a quintic 
nonlinear dispersion term to the parametrically driven 
complex Ginzburg-Landau equation suffices for obtain- 
ing a hysteretic nonequilibrium fsing-Bloch transition 
(HN1BT) , a phenomenon first observed in [j| . As a fifth- 
order nonlinearity represents the simplest, higher order 
correction to the usual complex Ginzburg-Landau, and 
the latter is of wide applicability in physical and chemi- 
cal systems, the HNIBT could be well observed in other 
systems. While we do not claim that this simple model, 
Eq. @ , represents an accurate description of the exper- 
imental system in || (the photorefractive nonlinearity is 
saturating) we note that it describes qualitatively the ex- 
perimental observations, even reproducing small details 
such as the shoulder in the field intensity in the case of 
the Bloch wall, Fig. (3), that has been repeatidly ob- 
served during the experiments in 
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Figure Captions 

Fig. 1. Homogeneous solution intensity as a function 
of detuning for 7 = 2, fi = a = — [3% = 1 and (3$ = 0.35 
(i), P 5 = 0.3928 (ii), and (3 5 = 0.45 (iii), which is an 
enlargement of a part of the positive detuning domain. 

Fig. 2. Homogeneous solution intensity for the same 
parameter values as Fig. 1 except ,05, which are marked 
in the figure. The different pattern domains are marked 
as IW (Ising walls), BW (Bloch walls), OW (oscillating 
walls), and P (patterns). The continuous (dashed) line 
indicates stable (unstable) homogeneous solution. 

Fig. 3. Intensity (full line) and phase (dashed line) 
spatial profiles of an Ising (a) and a Bloch (b) wall for 
the same parameters as Fig. 2(a) and the detuning values 
marked in the figure. 

Fig. 4. Intensity (full line) and phase (dashed line) 
profiles of an oscillating wall at three different instants of 
time (see text). The parameters are the same as in Fig. 
2(b) and v = 1.29. 

Fig. 5. Velocity of the domain walls as a funtion of 
detuning for the same parameter values as in Fig. 2. The 
gray areas mark the velocitie of the oscillating walls. The 
arrows mark the transition from Bloch walls to oscillating 
walls when a decreasing detuning scan is carried out. 

Fig. 6. Velocity of the domain walls as a function of 
detuning for 7 = 2, ji = a = 1, /?3 = — 1.5, and /?5 = 0.6. 
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